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Abstract 

Models of intermittent behaviour are usually formulated using a 
set of multiplicative random weights on a Cayley tree. However, inter- 
mittency in particle multiproduction from QCD jets is related to frag- 
mentation of an additive quantum number, e.g. energy-momentum. 
We exhibit the non-trivial stochastic mapping between these additive 
and multiplicative cascading processes. 

1. Introduction: Intermittency and its formal description in terms 
of mutiplicative cascading 

When Andrzej Bialas and I, 10 years ago, were puzzled by the mysteri- 
ous multiplicity fluctuations observed in an ultra-energetic cosmic ray event 
(JACEE collaboration), we were far from realizing that a systematic study 
of multiparticle production processes would emerge from considering the fac- 
torial moments of the multiplicities distributions. The proposal we madeS 
was guided only by two seemingly reasonable requirements: 

i) obtaining a non-subjective measure of fluctuations using the moments 
of the multiplicity distribution, 

ii) eliminating (as much as possible) the obvious source of statistical fluc- 
tuations due to the finite number of particles produced by event. 

We thus proposed the measurement of normalized factorial moments. 
Under the (strong!) hypothesis that the statistical noise was simple, i.e. 

*Written in honour of Andrzej Bialas for his 60th birthday and ever young passion for 
Physics. 
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Poissonian in the following formula (it was assumed Binomial in the case of 
the JACEE event because of the fixed multiplicity), one writes: 
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where m labels a (suitably chosen) piece of available phase-space for the 
reaction products, n is the number of particles registered in the bin for one 
event and p is by definition an associated continuous dynamical variable 
corresponding to the local multiplicity density. 

It is quite clear that the absence of fluctuations other than Poissonian 
would give F q (m) = 1. Moreover, varying the binsize m allows one to look 
deeper into short-distance fluctuations without being hidden by the set of 
increasingly erratic statistical fluctuations. In the same spirit, one would 
expect to observe dynamical objects, such as clusters, resonances, jets - by 
some signal in the range of scales governed by the size of these objects. 
However, Andrzej and I were soon confronted with a question of interpreta- 
tion: What means a continuous rise of the factorial moments with decreasing 
phase-space size? Indeed, no obvious scale emerged from the observed rise of 
factorial moments for the JACEE event (this feature was later on confirmed 
by systematic studies of various reactions). 

At that time having no experience of such a phenomenon in high-energy 
physics, we were inspired by the studies of turbulent hydrodynamic flowsSl, 
in building a mathematical toy-model (cf. a-model, in the terminology and 
the version of Ref.B which we used after modification o) reproducing a pat- 
tern of multiplicity fluctuations similar to what was surprisisingly suggested 
by the JACEE event. 

In its simplest version, for one event, the model is given by an (infinite) 
set of randomly independent numbers {W s } located along the branches of a 
Cayley tree structure (see Fig. la), where s (1 < s < v) is a branching step 
while k s (1 < k s < K) denotes which of the K branches at a vertex bears 
the specific weight W s . The multiplicity density profile is obtained for each 
individual bin by the following multiplication rule: 



where the path [m] = {ki, k s , ...k v } is uniquely associated to the phase- 
space bin m. Using the statistical independence of the weights and assuming 



p[m] = n W * 
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many steps v 3> 1, one gets: 



^4 = {W«V = (M)^, (3) 



where M (= X") is the total number of bins and, by definition: 

{W 9 } r = J r{W) W q dW (4) 

where r(W) is a normalized weight probability distribution satisfying {l} r = 
{iy} r = 1. As can be inferred from formulae (1-3), the multiplicity moments 
of an a-model (once the statistical noise is suitably de-convoluted ) show up 
whith a power-law dependence on the number M of phase-space bins, e.g. 
the resolution with wich one examines the system. This behaviour is charac- 
teristic of the phenomenon of intermittency in hydrodynamical turbulence; 
it appears as a consequence of the random-cascading multiplicative property 
of the a-model. 

Only a few years later, the ubiquous character of intermittent behaviour 
was recognised. In particular, an a-priori different type of cascading be- 
haviour, -though much more familiar to particle physicists-, was suspected^ 
and then shown® to possess intermittency properties. Interestingly enough, 
it is a direct consequence of Quantum Chromodynamics -the field theory of 
strong interactions - for the phase-space structure and development of quark 
and gluon jets. In this framework, the cascading mechanism can be called 
additive, since at each elementary vertex energy-momentum is conserved and 
fragmented among the decay products along the cascade. This is to be dis- 
tinguished from the a-model, for which local densities are multiplied during 
the cascading process and thus not additively conserved . It is the subject 
of the present paper to exhibit the transformation which asserts the equiv- 
alence between multiplicative and additive cascading mechanisms and thus 
their identical intermittency properties. 
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2. From a local to a global description of intermittent cascading 



As is explicited by the relation (3), the multiplicity density moments 
of the a-model unravel the structure of fluctuations in the local limit (i.e. 
short distance of order 1/M). Our aim is now to look for the system as a 
whole, i.e., its global description. For this sake it is convenient i' to introduce 
the (random) Partition function Pf(q) and an associated generating function 
Z v (u); One writes 

1 M / W \ q 

PM - 17 £ pU = £ n Cjf) ; Z M = (e- uP ^), (5) 

1V1 m m =l s =l v xv ' 

where one includes in the computation all the paths of the Cayley tree for v 
cascading steps, see Fig. la. Note that a thermodynamical formalism can be 
usefully introduced^ where W s /K acts like a Boltzmann weight and q as an 
inverse temperature. 

Interestingly enough, Z v {u) is known in Statistical Mechanics to obey a 
master equation^. Using an iterative procedure ( v — > u+1 ) adding one new 
step at the beginning of the cascade and also the statistical independence of 
distinct sub-branching processes (and after some work), one obtains: 

= (4)},- < 6 > 

The compact formula (6) will be at the root of the mathematical trans- 
formation between the additive and multiplicative versions of intermittent 
cascading. Let us for convenience introduce some interesting extensions of 
(6). 

Extension 1: Inserting Pf(q) instead of -P/(l) in the generating function 
(5), one gets a g-dependent master equation: 



Z v+ i(u,q) = 




which exhibits a scaling behaviour in terms of q (or temperature 0.) Note 
that the equation derived from (6-1) for the first moment reproduces the 
intermittency property (3). 
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Extension $ 2: The master equation can be easily extended to random- 
cascading processes including also random-branching (see Fig. lb). Introduc- 
ing a time variable t, an e probability of branching between t and t + e, a 
change of variables — > t ; + 1 — > t + e, and going to the limit e — > one 

gets 



Here too, the reader can check that the first moment equation obtained by 
functional derivation with respect to u leads to an intermittency property 
(3) (with a modified exponent). Note that Eq. (6-2) takes the familiar form 
of a gain-loss formula. 

Extension # 3: Without modifying the basic properties of the cascading 
process, one may introduce a generalized distribution r(t; Wi, Wjj, Wk)- 
Using a suitable change of variables, and after some transformations of func- 
tions and variables one writes the following master equation: 



We shall soon recognize that equation (6-3) provides a generic form of the 
additive cascading model of a jet based on QCD. 

3. Additive vs. Multiplicative cascading models 

In particle physics theory, however, intermittent behaviour has not been 
found directly under the form (6-1,3) of a multiplicative cascading process. 
It appears in the study of the multiplicity of gluons and quarks associated 
with an energetic jet in the framework of the resummed perturbative expan- 
sion of Quantum Chromodynamics. It has been for instance applied for the 
decays of Z°s into quarks and gluons. In the leading-logs approximation of 
perturbative QCD for jet calculations, one writes!! (for gluons): 




(6-2) 




(6-3) 



dZ 



(Q,u) = ^Jdz$(Q,z) [Z(Qz,u)Z(Q(l-z),u)-Z(Q,u)} (7) 
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where Z (Q,u) is the generating function of gluon multiplicity factorial mo- 
ments from an initial gluon jet characterized by the virtuality Q; $ (Q, z) is 
given in terms of the renormalized QCD coupling constant a s and of the triple 
gluon Altarelli-Parisi-Kernel Pq G (including quarks would transform (7) into 
a two- by-two matrix form without changing our main conclusions). Note that 
the important QCD property of angular ordering allows the derivation of the 
multiplicity distribution for any subjet of initial energy E and conical aper- 
ture G starting from the same function Z (Q = EQ,u) . As an important 
consequence, not only global but also local properties of multiplicity distri- 
butions of a QCD jet are determined by the solution of the equation (7). In 
particular angular ordering leads to the property of angular intermittency. @ 

As is explicit in formula (7), when compared to formulae(6-l,3), the gluon 
cascading process is generated by the fragmentation of energy-momentum 
between gluons at the vertex, using the energy-fraction variable z. Moreover, 
identifying the "time variable" with In Q, we find that Eq. (7) is not defined 
with equal-time observables (on contrary to Eqs. (6)), due to the mismatch 
between virtuality of a gluon and energy sharing. 

Yet, the two approaches are equivalent, as can be inferred!^ from a 
crucial property of the multiplicity distributions, namely KNO scaling @. 
This scaling property is verified in QCD@ within the same conditions as Eq. 
(7). One may write, at least at high enough virtuality^; the following scaling 
relation: 

Z(Q,u)=((u(n) q ), (8) 

where (^)q is the average multiplicity at virtuality Q. Let us insert the KNO 
scaling relation (8) into (7). We get 



dZ 
d\nQ 



^Jdz$(Q,z) x 



X 



C (« (™>q z ) C (u (n) Q(1 _ z) ) -((u (n) 



(9) 
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which, by a suitable change of variable can be cast into the following equiv- 
alent form: : 



1 Thcrc could appear some problems near the boundaries of the z-integcation in Eq.(7). 
However a check of validity can be made using Monte-Carlo simulations U . 



6 



dZ {Q,u) 
dlnQ 



J dW 1 dW 2 f{Q;W 1 ,W 2 ) x 

(10) 

x \z [Q, u^j Z (Q, u^j - Z (Q, «)} , 



by defining: 



Wl = 2^;W 2 = 2 (ii) 
Wq " Wq 

and 

r (Q; ^, W 2 ) = $ (Q; z) (W 2 - W 2 (Q; W,)) , (12) 

where z (Q; W\) and (Q; Wi) are given in terms of the functional form of 
(n)g . This functional form is obtained by solving the linear equation coming 
from Eq.(7) for the first moment, i.e. the first derivative with respect to 
u. The solution of this equation is in general much simpler than for the 
generating function or higher moments. Note that a solution of this simpler 
equation is sufficient to define the transformation we look for (and thus to 
derive the intermittency properties). 

Equations (10-12) are to be compared with the generic equation (6-3). 
More precisely, the change of function and variable 

J JdW 1 dW 2 f(Q;W 1 ,W 2 f) 1 f(Q;W u W 2 ) -> r(Q;W u W 2 ) 

(13) 

J JdW 1 dW 2 f(Q;W 1 ,W 2 )^ dlnQ -> dt, 

gives the equivalence of gluon cascading with a random branching, random 
cascading, multiplicative process. It amounts to choosing in Eq. (6-3) K = 2 
and the specific probability distribution r, see (13). In a sense, the multiplica- 
tive formulation provides a statistical description of the multiplicity density 
distribution at equal time (or virtuality, or jet aperture), while the additive 
one describes the sharing of energy-momentum. Both descriptions are inti- 
mately connected by the KNO scaling property. It is worth noticing that this 
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specific stochastic process is, strictly speaking, of semi-random type since the 
value taken by the random weight W\ fixes the other one W 2 at the same 
vertex. 



4. Application to QCD gluon cascading 



Let us consider QCD cascading at the leading logarithmic approximation, 
by keeping into account the effect of energy-momentum conservation at the 
vertex. This problem has been recently raised for the fluctuation pattern and 
solved^. In terms of global observables, it amounts to solving equation (7) 
with the kernel $ (Q 2 , z) = 2^j z where 7q = 2a s N c /iT is kept fixed for sake 
of simplicity. Using the formulation of Eqs. (10-13), one gets: 



dZ (Q,u) 
d\nQ 



dW 1 

• i x 
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(14) 
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where the relations: 
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are coming from the expression of the mean multiplicity (^)q oc Q 70 ; by 



inversion, one finds: 



'Wi \ 1/7 ° 

z = z (Wi) — ( 



Wo 
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Wi \ 1/7 °] 70 

~2 



r(Q;W l ,W 2 )= lQ — 8\W 2 -2 



~2 



(16) 



A few remarks are in order about the master equation (14). 

i) Comparing (14) with Eqn.(6-3) shows that Z is the generating function 
for a random-branching process with triple vertex 1 — > 2, (up to transfor- 
mations like (13)). At each vertex, one branch corresponds to a randomly 
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chosen weight W\ with the probability law 70/W1, see expressions (16). The 
second-branch weig ht is then determined by W 2 = 2 [l - (Wi/2) 1/7o F° . In- 
terestingly enough, the local constraint of energy-momentum conservation 
along the additive cascading process is transposed into a specific stochastic 
law at the multiplicative vertex of the density cascading process. 

ii) In the approximation W\ < 2 ; 70 C 1, then (^f 1 ) ° "C 1, and one 
can write a simplified version of equation (13) namely: 



<91nQ IU J Wi V V 2 

Equation (17) has been found® equivalent to the double-leading-log approx- 
imation of QCD for the jet process, where one can neglect the recoil effect 
upon the leading parton (quark or gluon). However, sizeable values of 70 
(which is still of order .5 at LEP energies) do lead to substantial modifi- 
cations of the multiplicity distributions'!!-} and their fluctuations in phase- 
space0 due to energy-momentum conservation. 
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Figure Caption 

a) Fixed branching 

Each path [m] corresponds to a bin in phase space. The fluctuation density 
in the bin is described by the product of random weights W along the path 
[m]. In the case of "semi-randomness" the decay is assymetric at each step: 
one weight W\ follows a random law, while W% is a function of W\ for the 
same vertex (W\ or W 2 can be randomly left or right). The tree can be 
separated into 2 branches. Each of these branches defines in average Z v of 
the Partition function, the whole tree itself corresponding to Z u +i- 
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b) Random branching 

Same structure as in Figure l.a, but with small increments e of the num- 
ber of steps v and random-branching with uniform probability e. 
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